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We study the correlation between spin-independent and spin-dependent scattering in the context 
of MSSM neutralino dark matter for both thermal and non-thermal histories. We explore the 
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Q-i generality of this relationship with reference to other models. We discuss why either fine-tuning or 

r^, numerical coincidences are necessary for the correlation to break down. We derive upper bounds 

^ on spin-dependent scattering mediated by a boson. 
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I. INTRODUCTION 



There is compelling evidence for dark matter (DM) from a variety of astrophysical 
observations: rotation rates of galaxies, gravitational lensing, and temperature perturbations 
of the cosmic microwave background (CMB). The WMAP experiment, which measures the 
CMB, provides the best value for the DM relic density [Ij, 

^]j3M/i^ = 0.1131 ±0.0034. (1) 

While there are potential hints of DM in cosmic ray anomalies {e.g. PAMELA [2] and 
Fermi [3]) and from the annual modulation signal seen by the DAMA/LIBRA collaboration 
[U |5], there is currently no measurement which gives conclusive information about the DM's 
properties. The identity of the particle responsible for I^dm remains a mystery. 

One DM candidate worthy of special attention is the Weakly Interacting Massive Particle 
(WIMP). The strongest motivation is the so called "WIMP miracle" - a thermally produced, 
stable particle with a weak scale mass and perturbative couplings will freeze-out with a relic 
density of the right order of magnitude to constitute the DM [6j. When added to the 
expectation of new physics at 0{TeV) responsible for the stabilization of the weak scale, 
the WIMP paradigm becomes even more compelling: the new physics often includes novel 
symmetries, which result in at least one particle being stable on cosmological timescales. If 
this particle is weakly interacting then it can be WIMP DM. 

However, this picture turns out to be an over-simplification. Not just any weak scale stable 
particle will do. If the DM is weakly interacting in the strictest sense - i. e. has full-strength 
SU{2) L^U (l)y gauge interactions - then DM may be excluded by existing direct detection 
(DD) experiments. In particular, a weak-scale Dirac (vector-like) fermion, xd? with SU{2) 
interactions (which encompasses the simplest DM model of all, a Dirac neutrino), feels the 
weak force via the operator: 

aector = (XD7^XD)^2. (2) 

When the coefficient of this operator is typical in size, namely 0{g/ cos^^u), where g is the 
SU (2) coupling constant and 9^ is the weak mixing angle, it leads to a huge DD signal 
- experiments constrain the DM mass to be greater than 50 TeV [7]. Furthermore, the 
thermal relic density for a 50 TeV Dirac neutrino will be far too large to explain the WMAP 
measurement. Thus, DM at the weak scale requires a strong suppression of this operator. 
In fact, it is straightforward to eliminate it entirely. If x is a Majorana spinor, the operator 
ix 7^ x) identically vanishes due to the properties of Majorana bilinears. The DM may be 
Majorana if an SU{2) singlet Majorana fermion mixes with a Dirac state. This mixing can 
only be accomplished via SU{2) breaking in the WIMP sector, i.e. through a Higgs boson 
vacuum expectation value (vev). Then the resultant DM particle has a non-zero coupling to 
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a Higgs boson, h, and the dominant scattering process is due to the following operators: 




(3) 
(4) 



In a multi-Higgs boson theory, h need not be the Higgs boson of the Standard Model (SM), 
but even in these theories, there often is a Higgs boson that has SM-like properties. We will 
explore the impact of these operators on Spin-Independent (SI) and Spin-Dependent (SD) 
scattering off of nuclei, paying particular attention to the expected correlation between the 
rates at these two types of experiments. 

While we perform most of our analysis in the context of the Minimal Supersymmetric 
Standard Model (MSSM) (for a review of the MSSM, see [8]), we reference other models 
where appropriate to emphasize the generality of our arguments. We will review the assertion 
that post-LEP (largely due to the constraints on the chargino and slepton masses), one may 
consider a mixed or "well-tempered" neutralino as a likely DM candidate, if it is thermally 
produced [9]. We will show that in this case, light Higgs boson and exchange will 
generically lead to a signal in the next generation of SI and SD experiments. 

A thermal history for the WIMP is not the only possibility. For example, non-thermal 
mechanisms may populate the DM {e.g. through the decay of a modulus or gravitino [TO]), 
or the DM can be overabundant and subsequently diluted by extra sources of entropy. These 
options allow a WIMP with a wider range of properties, since the annihilation rate is not 
fixed by the thermal history. In what follows, we do not rescale DD signals to the (too- 
low/too-high) thermal relic density. In all cases, we assume that the WIMP constitutes the 
total DM density, determined from astrophysical measurements to be pdm ~ 0.3GeV/cm^. 
We will be clear when we are making the assumption of a thermal history. For the purposes 
of this study, a "thermal" WIMP is one whose thermal relic density is within the generous 
range ±3cr of the WMAP measurement given in Eq. ([T|. We also note that more recent 
determinations favor a slightly larger value: p^u = 0.39 — 0.43 GeV/cm^ [111 [12] • This 
would extend the reach of the direct detection experiments by a factor of ~ 4/3 and probe 
more parameter space. An accurate determination of the local DM density is important for 
an accurate measurement of the DM DD cross section. 

Related results already exist in the literature, including some comprehensive numerical scans. 
However, we find that often the (simple) underlying physics is left obscure. We hope to make 
clear the expected size of various contributions to DD and the relationship to the assumption 
of a thermal relic abundance. Assuming there are no conspiratorial cancellations, these 
typical sizes represent important targets for DD experiments. 

There is an overwhelming literature in existence on the subject of DD, see reviews [TSj [H] 
and references therein. Of particular interest to us is the relationship between the size of the 
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SD and SI signals, which has recently been explored in |15l4T8] . 

In the next section, we begin by discussing the current experimental status and then make 
naive estimates for the SI and SD DD cross sections from h and exchange respectively. In 
Sec. in we lay out the specific structure of the SI and SD operators in the MSSM and estimate 
the naive size of the SI and SD cross sections. Then in Sec. IV we review the argument for a 
well-tempered neutralino and discuss some alternatives. Sec. |V] concentrates on illuminating 
the expected size of the SD cross section for mixed DM models with various restrictions. In 
Sec. |VT]we describe the conditions under which SI and SD signals in the MSSM are expected 
to be correlated. Technical results are relegated to three appendices. 



II. DIRECT DETECTION PRELIMINARIES 



The interactions in Eqs. (|3]) and (|4]) lead to SI and SD elastic signals in DD experiments, 
respectively. In Fig. [T]we have plotted the current experimental limits for SI and SD DD. 
Currently, the state of the art SI experiments are CDMS [H] and XENON [201 ■ XENON 
constrains a^f < 4.5 x 10^® pb for = 30 GeV. After combining their most recent run 
with previous data, CDMS-II has a 90% CL bound of 3.8 x 10"^ pb for a WIMP with a 
mass of 70 GeV In the most recent data set, two tantalizing events were seen, but it 
is premature to attribute these to signal. In any case, XENONIOO expects to place a limit 
on the order of a^f ~ few xlO~^ pb by early 2010. Thus, we will consider SI cross sections 
greater than 5 x 10~^ pb as potentially probeable in the short-term, and hence "large." 

There are two ways the SD cross sections are constrained. The first is via DD experiments. 
The current best bound on the SD DM-proton interaction comes from the KIMS experiment 
[22], cr^D < 1-6 X 10"^ pb for = 70 GeV; the best bound on the SD DM-neutron 
interaction coming from the XENON experiment, a^^ < 6 x 10"'^ pb for = 20 GeV, with 
the strongest bounds for masses of (9(10) GeV coming from PICASSO [22] ■ There are also 
bounds from DM capture in the sun, assuming (as is the case in the MSSM) that the DM has 
annihilation products which give rise to relatively hard neutrinos. Assuming annihilation of 
the DM to bosons is appreciable (as is appropriate for much of the parameter space 



considered here, see Sec. IV), IceCube [24] places very strong bounds for masses above 250 



GeV with the strongest bounds coming at 250 GeV, cxg]^ < 3 x 10""^ pb. At present, no 
limits exist from IceCube below this mass. For smaller masses, the best limits of this type 
come from SuperK [23], < 10~^ pb above > 20 GeV. 

Perhaps within the next two years [26j, the COUPP [27] and PICASSO [2H] experiments 
will take data with a projected sensitivity to SD scattering of cTg^ ~ 10~^ pb. They will also 
have sensitivity down to much lower masses than the neutrino experiments. The XENON 
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data will probe a^^ ^ A x IQ-^ pb for a 30 GeV WIMP. A 1 ton COUPP-like proposed 
experiment [IT], might ultimately probe values as low as 10^^ pb. The DeepCore extension 
to the IceCube detector should be able to extend down to the 10~^ pb level with 5 years 
of data [29]. Bounds from neutrino experiments can be avoided if particular final states 
dominate WIMP annihilation, e.g. generation quarks, though this does not happen in 
the MSSM. We consider SD cross sections greater than 10""^ pb as potentially achievable in 
the short-term, and hence "large." 
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FIG. 1: Current bounds on SI (left) and SD (right) DM-nucleon cross sections. The COUPP and 
XENONIOO projected SD bounds are only estimates - we have scaled the current exclusion curve 
of COUPP by a factor of 10~^ [30] and the current SD exclusion curve of XENONIO by the factor 
which scales the XENONIO SI limit to the XENONIOO SI limit. 
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A. Spin Independent 



The operator responsible for SI DM-nucleus interactions is 

of = c,{xx){qq), (5) 

where x is the DM and g is a quark. Taking the expectation value of this operator between 
two nucleon states (A^ = p (proton) or n (neutron)) determines the effective interaction of 
the DM with a nucleon, 

{N\m,qq\N)=mj,f!fl\ (6) 

where the nuclear matrix element f^^^ is determined in chiral perturbation theory from the 
pion nucleon-scattering sigma term. The coefficient of the effective DM-nucleon interaction, 
fN{,XX){.N N), is given by 

niN ^-^ rrin 27 m„ ^ 

q=u,a,s ^ q=c,b,t ^ 

where = 1 — X]g=u d s frq ^"^^ ^ ^^^^ refers to Higgs boson exchange [31] . 

The nucleon-Higgs interaction is coherent over the nucleus [32] resulting in the well known 
enhancement for SI cross sections. To compare between experiments using different nuclei, 
the elastic scattering cross section is normalized to a per nucleon value [13] : 

asi(x iV ^ X iV) = ^ ml j^{Zf, + {A- Z) ^)^ (8) 
where is the reduced mass between the DM and a nucleon. 

We use the DarkSUSY package for numerical analysis [33], so for analytic estimates we will 
use the same values for the nuclear matrix elements, namely 

f!^l = 0.023 = 0.034 ff^] = 0.14 f^^^ = 0.803 
= 0.019 f^""] = 0.041 = 0.14 = 0.800. 

Since these values are derived from the pion nucleon scattering sigma term, their error bars 
are correlated. 

If the Higgs boson, h, that mediates the interaction between the DM and the nucleon is 
SM-like, the coefficients Cq are given by 

1 
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where Uq (y^) is the Yukawa couphng for the quark (DM) and ruh is the Higgs mass. The 
per nucleon cross section is then 

MxN^ XN)^5X 10-« pb (^)' (^Hi^y (SI typical). (10) 

Estimates based on recent lattice simulations seem to favor smaller values for the nuclear 
matrix elements [3l]. If these lattice results are correct, the dominant contribution to the 
SI scattering cross section would be due to the heavy quark content of the nucleon (since 



fxG 1 in the limit of small fq^^) and the coefficient in Eq. (10) would be replaced by 
2 X 10^*^ pb. In cases where q ^ Cu, which can occur in models with multiple Higgs 
bosons such as the MSSM, then uncertainties in the f!j^^ can lead to as much as an order of 
magnitude variation in asi{x^ X^) [35] . 



B. Spin Dependent 

The operator responsible for SD DM-nucleus interactions is 

O^"" = dq{xr^'x){qi,l'q)- (11) 

Taking the expectation value of this operator between two nucleon states allows us to find 
the effective SD interaction of the DM with a nucleon (A^ = p (proton) or n (neutron)), 

(iV|g7,75g|iV) = 2.WAg(^), (12) 

where s^/^^ is the spin of the nucleon and the Ag*^^-* are extracted from polarized deep elastic 
scattering. The coefficient of the effective DM-nucleon interaction, 2aN{x 7'^7^ x) (^ ■Sm^'' N) , 
is given by 

aN= Yl dqAq^'^K (13) 

q=u,d,s 

The elastic scattering cross section quoted by the experiments is between the DM and a 
nucleon which is given by 

asBixN ^xN) = -mla%, (14) 

IT 

where rrir is the reduced mass between the DM and a nucleon. 

Again we follow DarkSUSY and use the following values for the SD calculations, 

a2') = 0.77 aJ') = -0.40 Ai^) = -0.12 

aW = -0.40 A;;"^=0.77 Ai^) = -0.12. (15) 
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The prediction for SD scattering is somewhat more robust to variation in the hadronic matrix 
elements than the SI case: the uncertainties in these values can lead to (9(30%) variation in 
the SD cross section [35] . 

If the SD interaction is mediated by the boson, then the coefficients dq are given by 

where Qz-dm parametrizes the coupling of the DM to the and Cw = cos 6*^. For 
concreteness, (and since it is relevant for calculations of solar capture) when we quote values 
for SD scattering we will focus on the cross section off of protons. For SD scattering mediated 
by the Z^, the neutron scattering is 0(20%) smaller. The SD cross section is 

c^sd(xP ^ Xp) ^ 4 X lO-^b (^ff^) (SD typical). (17) 
In the next section, we discuss the form that Qz-dm takes in the MSSM. 



III. DIRECT DETECTION OF NEUTRALINO DARK MATTER 



The best DM candidate in the MSSM is the lightest neutralino, which is an admixture of 
Bino (B), Wino (W), and the up and down-type Higgsinos {Hu and Hd). The stability of the 
lightest superpartner (LSP) is guaranteed by i?-parity, which is introduced to avoid proton 
decay. The neutralino mass matrix is given by 



M 



f Ml 
M2 

-mz 0/3 mz Cyj cp 

\ mzSujSis -nizCwSp 



-mz cp mz 
mz Cyj Cp -mz Sp 
-/i 







where Mi is the Bino mass, M2 is the Wino mass, /i is the Supersymmetric (SUSY) Higgs 
boson mass parameter, mz is the Z^ mass, {3 = arctan(t>„/t>rf), are the up and down-type 
Higgs boson vevs, = sin6'^, = cosOw, = sin/3, and = cos/3. 

The composition of the lightest neutralino, which we denote x, is specified by 

X = ZbB + ZwW +ZH,Hd + Zh^ Hu. (18) 
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If squarks are heavy, the only potentially sizable contributions to SI DD are from both CP- 
even Higgs bosons, h and H, where m/j < mn- We comment on the typically subdominant 
squark exchange contributions in Appendix |Aj The Higgs boson exchange contributions are 
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where Cu4 are the SI operator coefficients given in Eq. (|5]), is the SU (2) gauge coupling, mw 
is the mass, = tan^^, a is the Higgs mixing angle, Cq, = cos a and Sa = sin a. In the 
decoupling {rriH oo and a — > 7r/2 + /3) and large limits, these expressions simplify: 
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(19) 
>(20) 
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(21) 
(22) 



where we have only kept the tj^ enhanced contribution from H . We will use these expressions 
below in Sec. VI when analyzing the allowed suppression of the SI cross section. 

The lack of an observation of a Higgs boson at LEP makes it likely that we live in at 
least a moderate tp regime (so that the tree-level contribution to the Higgs boson mass 
irih = mz cos 2/3 is maximized), and constraints on the mass of the charged Higgs from flavor 
experiments point to the decoupling limit. Therefore, Eqs. (21) and (22) are particularly 
useful for estimating the expected size of scattering. In Sec. IV we will argue for the typical 
size of the various neutralino mixing angles which lead to SI cross sections of the order. 



cr, 



MSSM 
SI 



5 X 10"^ pb 



115 GeV\' 



JW 



tyj Z 



B) ^Hu 



0.1 



(MSSM: SI typical), (23) 



where we have used Eqs. (21) and (22) and taken niH ^ oo. 



In the heavy squark limit, contributions to SD DD come from Z^ exchange. Since the Bino 
and Wino are both SU{2) singlets, they do not couple to the Z^. Therefore, SD is controlled 
by the Higgsino content of the WIMP. The exchange contribution takes the form: 

,2 



d„ 
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- \Z 



(24) 



A non-zero Higgsino component (so that Zh^ ^ 7^ 0) is insufficient to ensure a non-zero SD 



coupling. If Ml, M2 — )■ 00, so that a pure Higgsino is recovered, \ZhJ 



1/^2, and 
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the SD coupling vanishes. Instead, the Higgsino forms a Dirac state, and the large vector 
scattering of the Dirac neutrino is recovered. Hence, mixing with B and/or W (so that 
\ZhJ 7^ l^Hdl) is required in order for the dq^s to be non-zero. This requirement also implies 
a non-zero SI cross section, giving the correlation demonstrated below. 



The typical cross section for SD DD in the MSSM (again see Sec. IV) is given by 



<D^^(XP ^ XP) ~ 4 X 10"^ pb ( ) (MSSM: SD typical). (25) 





2 _ 






0.1 J 



There are reasons to expect the squarks do not make a sizable contribution to the DD cross 
sections. In the MSSM, satisfying the LEP bound on the Higgs boson mass requires large 
radiative corrections from the stop loops. This implies that at least one stop must have a 
TeV scale mass. Renormalization group flow tends to make the third generation sparticles 
lighter than the partners for the first and second generations. Therefore, it is plausible that 
squark contributions to DD scattering are negligible since only the first and second generation 
squarks contribute (see Appendix|A]for details about squark exchange). For concreteness, in 



all scans below we take the scalar superparters to be (9(2 TeV). This is also why Eqs. (23) 



and (25) are expected to be good approximations. For a study which focuses on the effects 



of light squarks, see 



IV. THE ARGUMENT FOR A WELL-TEMPERED NEUTRALINO 



Arkani-Hamed, Delgado and Giudice P] argued that when one takes the LEP limits on 
charginos and sleptons into account, a pure neutralino {i.e. composed of only one gaugino 
eigenstate, usually taken to be Bino) is no longer the "natural" MSSM DM candidate, at 
least when one imposes the requirement of a thermal cosmology. They claim that one should 
instead consider a mixed neutralino, which they have dubbed "well-tempered." Since the 
relic density of mixed DM is set by annihilations to W~ (and 1 1 when kinematically 
allowed) there is a further condition that > my/- Hence, we will impose this requirement 
when we refer to "thermal" DM in the analysis that follows. In what follows, we review their 
argument and then discuss some non-thermal options. Note that SI DD has previously been 
studied for well-tempered models [SHlEn], but no dedicated SD study exists. 
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A. Thermal history 



We begin by considering the thermal history of a nearly pure Bino. If one does not allow 
for CO- |^0] - H2] or resonant |10l H31 - H6] annihilations, then Bino freeze-out is controlled by t- 
channel sfermion exchange. One can show [9] that in order to produce the observed DM relic 
density, the sfermion must be ^ 110 GeV. Since the LEP limits on sfermions are (9(100 GeV), 
there is only a small experimentally allowed window for thermal Bino DM. 

Either co-annihilations {e.g. with the stau or stop) or resonant annihilation through the 
pseudo-scalar Higgs {A^) also allow dominantly Bino DM. However, both of these options 
involve numerical coincidences. In the first case the Boltzmann factor will exponentially 
suppress the density of the would-be co-annihilator unless exp(— AM/T/) is 0{1), where 
AM = m-NLSP — ''^x' ''^NLSP is the mass of the next-to-hghtest superpartner, and T/ is the 
DM freeze-out temperature. Since Ty ^ m^/20, this requires a mass degeneracy, AM, of a 
few percent. To realize the second case requires a precise relationship between and rriA- 
When < mw, the Z° or h poles may be used to achieve the correct relic density, which 
requires a similar numerical conspiracy. 

Located at the other extreme, far away from the pure Bino, is a pure Wino or a pure Higgsino. 
In these cases, the requirement of a thermal relic abundance fixes the mass to be (9(2.5 TeV) 
and (9(1 TeV) respectively. Thus, to realize either of these cases implies /i > (9(100 GeV). 
Since, in the MSSM, the mass is given by 

^ = -H' + ^^^^pf^, (26) 

where ^ are the Higgs soft-mass squared parameters, this requires a substantial fine- 
tuning between /x^ and irn^ ^ in order to reproduce the measured mass of 91 GeV. 
Therefore, the desire to alleviate fine-tuning in this expression leads to the requirement that 
fi ~ (9(100 GeV). This will also naively lead to well-tempering since the neutralino mixing 
is proportional to rriz / fi- Though the accuracy of the current measurement of the DM relic 
density (see Eq. ([T])) requires a precisely determined neutralino composition, one can easily 
reproduce the DM abundance for any mass of (9(100 GeV). The Bino/Higgsino mixed LSP 
as a good thermal WIMP was pointed out in studies of the focus point region on the MSSM 

miiiB]. 

A Higgs boson mass above the LEP bound requires large radiative corrections from a stop 
squark. This implies that the scale for these particles, msusY, should be around a TeV. These 
states yield additive corrections to m^j^^, proportional to mgusv Hence, even in the case 
when /i ~ (9(100 GeV), there will naively be fine-tuning between these corrections and the 
bare value of w?jj in order to reproduce mz- Solutions to this "little hierarchy problem" 
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have been proposed within the MSSM {e.g. [19]) - we will ignore this type of fine-tuning in 
our arguments, focusing instead on the model independent tuning explicit in Eq. (26). 



B. Non-thermal options 

A thermal history is not the only way to achieve the correct DM relic abundance [TU] ■ It has 
even been argued [SU] that there is a "non-thermal WIMP miracle" when there exist TeV 
scale states which decay to the DM via Planck suppressed operators. For example, a heavy 
gravitino (or string-theory moduli fields) can live long enough to dominate the energy density 
of the universe. Then when these states decay, they will produce superpartners which will 
decay down to the lightest neutralino, resulting in a neutralino relic density. This relaxes 
the relationship between the mass/composition and relic density of a neutralino. 

A variety of other options have been proposed. Models where the energy density of the 
universe at the epoch of DM freeze-out was dominated by something other than radiation 
were studied in pT] . Alternately, if the DM interacts so feebly that it never achieves thermal 
equilibrium, one can achieve the correct value of the relic density via "freeze-in" production 
[52] . Since the total energy density of DM is close to that of the baryons, one can construct 
models where the DM relic density is set by an asymmetry which is determined by the baryon 
asymmetry [S3]. In [M], it was shown that by varying the reheat temperature and allowing 
for non-thermal sources, any neutralino composition can result in the correct relic density. 
In [SS] , a low temperature phase transition in the early universe changes the DM properties 
after freeze-out. All of these options involve either non-trivial cosmological histories or other 
model building challenges. We will focus on the thermal - and hence well-tempered - case, 
with discussions of the deviations that arise when the thermal assumption is relaxed. 



V. SPIN DEPENDENT CROSS SECTIONS FOR MIXED DARK MATTER 

In the MSSM, the neutralino mass mixing can often be approximately understood in terms 
of a two state system: a Dirac Higgsino mixing with either a Bino or a Wino. Thus, to 
understand the physics of SD scattering via Z° exchange, it is useful to consider the simple 
"Singlet-Doublet Model" (SDM) for DM, where the singlet has the same quantum numbers 
as either a Bino or a Wino, and the doublets have the same quantum numbers as the 
Higgsinos: 

CsDM^ I^DDD + XhSD + X'h*SD + ^S\ (27) 
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Here D and D are a vector-like pair of SU{2) doublet fermions, 5* is an SU(2) singlet, h is 
the SM Higgs doublet, A (A') is the Yukawa coupling which leads to the mixing between the 
D (D) and S", /id is the vector-like mass for the D and D, and fis is the Majorana mass for 
S. For the purposes of SD scattering it is sufficient to replace h by its vev, {h) = v = 174 
GeV. The exchange of the uneaten component of h leads to SI DD. 

In the case where S plays the role of the Bino, the values of A and A' are constrained by 
the supersymmetric relations to be A f = —mz Sw Cj^ and \' v = —mz Sw Sj3, while in the case 
where S is the Wino, the values of A and A' are constrained by the supersymmetric relations 
to be Af = mz Cw Cjs and \' v = mz Cw Sjs. 

We now use this model to discuss the coupling of the boson to the DM in the MSSM. 
In Appendix |B] we discuss the diagonalization of the 3x3 mixing matrix of the SDM. With 
appropriate substitutions, these expressions correspond to either Bino/Higgsino (M2 — )• 00) 
and Wino/Higgsino (Mi — )■ 00) neutralinos. In these limits we can write down approximate 
expressions for the effective coupling of the DM to the Z^. When there are no degeneracies 
between parameters in the neutralino mass matrix and mz may be treated as a perturbation, 
we have (see p] and Appendix [B| : 



\ZHf-\ZH\'=l .-M^ -^'--I'l^l'l^l 



|2 



for |Mi 


1. I/^I. 




- Ml > mz, M2 -> 00 


for IM2 






- M2 > mz, Ml 00. 


do not 


occur 


in 


this limit. Instead, they are 



The largest values of \ZhJ — \ZhJ\ 
when two parameters of the neutralino mass matrix are degenerate. The reason is simple: a 
degeneracy allows a large gaugino-Higgsino mixing in spite of the relative smallness of the 
off-diagonal entries of the neutralino mass matrix (proportional to mz)- It should be said 
that there is no particular reason to believe that a precise degeneracy should occur, since 
and the gaugino masses are SUSY preserving and breaking respectively. However, since this 
case maximizes the possible signal at SD experiments, it is worth noting. In the presence of 
these degeneracies, we have (see Appendix [bJ: 

I — ^TlH — + 8^^^ I - 1^1 > T^z, Ml 00. 

Perturbing away from the limit of exact degeneracy gives corrections to these expressions of 
0{{Mi — Note that DM with a mixed Wino/Higgsino has a SD DD rate enhanced 

relative to a Bino/Higgsino admixture by the appropriate power of c^jsw = 1-8. 

What is the largest obtainable SD cross section in the MSSM? A numerical scan yields 

iZnf-lZnJ^ < 0.4^ (30) 
(f^SD^^) < 6 X lO^^pb (General MSSM, Non- thermal DM), (31) 
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when the squarks are heavy. This upper bound is largely a consequence of the LEP bounds 



on the chargino masses which force the mixing ~ to be less than one. Eq. (29 ) provides 

a good analytic understanding of this number - it comes within approximately 10% of this 
value. The deviation is due to mixing effects that occur away from the large Mi limit. 

In many models of SUSY breaking the relation Mi/ai = M2/a2 = M^/a^ holds. We refer 
to this condition as unified gaugino masses. Because this is equivalent to M2 ~ 2Mi at the 
weak scale, the LSP is mostly Bino and Higgsino. In this case, 

\ZH,\^-\ZHf < 0.32^ (32) 
(^SD^^) < 4 X 10"^ pb (Unified Gaugino Masses, Non - thermal DM). (33) 

Finally, for > mw-, a thermal relic density within ±3 a of the WMAP measurement 
implies an upper limit on the amount of Higgsino in the DM particle. Therefore, 

\ZHf-\ZHS < 0.24^ (34) 
(a|g^^)thermai < 2 X lO-^pb (General MSSM, Thermal DM) . (35) 

This result holds for the case with unified gaugino masses as well. Note that Eqs. ([3T|), (p3|. 



and (35) all occur for a DM mass of C(80GeV). 



To saturate the above bound {i.e. maximize ctsd for thermal, well-tempered DM) requires a 
Bino/Higgsino mixture (recall that dq vanishes for a pure Higgsino), with a negligible Wino 
contribution. The largest values of SD DD occur when the DM has the largest Bino/Higgsino 
mixing which happens for the lowest values of the DM mass. As the mass of the DM increases, 
a larger component of Higgsino or Wino is needed for the DM to efficiently annihilate down 
to the correct relic density, which in turn typically leads to a decrease in ctsd- 

As shown in Fig. [2| there is a tight correlation between the SD cross section and the DM 
mass, in the decoupling limit when there is gaugino mass unification and a thermal relic 
abundance. 

For low masses, the neutralino is well-tempered for low masses and as — )■ (9(1 TeV) the 
neutralino approaches a pure Higgsino. Examining Fig. [2| except for when the annihilation 
channel XX ~^ tt opens, ctsd is a smooth, monotonically decreasing curve. An experiment 
sensitive to cross sections of (9(10^'^ pb) will probe < 200 GeV. There is a spread in the 
points in this figure from the liberal range taken on the relic density constraint. For masses 
approaching (9(1 TeV), there is additional extent from the variation in the Bino content of 
the neutralino and from contributions from squark exchange. For masses at 1 TeV, ctsd goes 
from 10~^pb — )■ for Mi from 1300 GeV — )• 00. Note that the projected reach of a 1 ton 
version of COUPP is O{10~^) pb for = 1 TeV [l7j, which would probe the entire range 
of SD cross sections for neutralinos excepting a nearly pure TeV Higgsino. 
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Note that the imposition of the unified gaugino mass condition essentially imposes the 
requirement that there is a tiny Wino content in the LSP. The hatched region in Fig. |2] 
is filled in when non-unified gaugino masses are allowed. In this case, a thermal relic DM 
candidate can be obtained for a Bino tempered with Wino if Mi ^ M2, which implies that 
the SD cross section decreases, effectively filling in the region beneath the curve in Fig. [2j 
Note that when ctsd ~ C(10~^pb), there is additional model dependence since the squark 
contribution becomes important (see Appendix IaI). 



Thermal Relic Density 
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FIG. function of for points satisfying the relic density constraint. We have imposed 

gaugino mass unification and taken the decoupling limit. The shaded region above the dotted line 
corresponds to "large" SD and will be probed in the near term. The solid red line is the current 
bound from IceCube, assuming annihilation to W~ . The blue hatched region is filled in if the 
assumption of gaugino mass unification is relaxed. The sfermion masses are taken to be 0(2 TeV). 

Finally, we note that there is a region of well-mixed Higgsino-Wino near 2 TeV with a 



thermal abundance (where M2 ~ /i). In this case, the second line of Eq. (29) applies, and we 
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find an approximate SD cross section of 6 x 10~^ pb, perliaps able to be probed at a future 
1 ton COUPP-like experiment. Tliese are the neutralinos which account for the hatched 
region above the points in Fig. |2} 

Not only is the SDM a simplified system useful for understanding the physics of SD scattering 
in the MSSM, it is potentially of independent interest. The DM may be unrelated to the 



solution to the hierarchy and simply given by the Lagrangian of Eq. (27) [561 EZ]- Then the 



DD story is essentially unchanged except there is greater parametric freedom. 

For example, the Higgs boson mass is no longer fixed by SUSY. Then the only constraint is 
nih < (9(TeV) to unitarize scattering. For m/j ~ TeV, the SI DD cross section is at most 
10~^^pb which would not lead to a signal in the next round of SI experiments. While such 
a large Higgs boson mass is in tension with precision electroweak measurements, it could be 
reconciled with a contribution to the T parameter [58] in a way that factorizes from the DM 
phenomenology. 

If one allows for a non-thermal history, the freedom of the SDM allows off-diagonal 
parameters of the mixing matrix that give {ZhJ"^ — \Zhu\^ = 1- This maximizes the SD 
DD signal from exchange (o"g§^ 4 x 10^^ pb). Thus, the SDM with a non-thermal 
history predicts scattering anywhere up to (or even above) the current bounds. Requiring a 
thermal history limits the amount of doublet allowed in m^, decreasing crg§'^. For if a very 
large doublet component is chosen (in an attempt to maximize the SD cross section), the 
requirement of reproducing the relic density requires ns to be C(TeV). 



VI. SPIN INDEPENDENT VERSUS SPIN DEPENDENT 



When a Majorana fermion couples to the there is necessarily an interaction with a 
Higgs boson, which leads to SD and SI elastic scattering respectively. In the last section, 
we concentrated on the physics behind the size of the SD cross section. We now ask the 
following questions: what is the expected correlation between the SI and SD signals? Is it 
possible to make one large while the other nearly vanishes? 

Since rrih and mz are known in the MSSM, there exists a correlation between the SI and 
SD signals, at least in the limit of heavy sfermions and Higgs boson decoupling. For this 



region of MSSM parameter space, the SI and SD DD cross sections are given by Eqs. (23) 



and (25), where only mixing factors and the Higgs boson mass are left unspecified. The light 
Higgs boson mass is constrained to lie in the tight range 114 GeV < uih < 130 GeV, where 
the lower bound is due to the LEP limit and the upper bound comes from considerations 
of fine-tuning. For the SplitSUSY model - where the decoupling and heavy sfermion limits 
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certainly apply - the Higgs boson mass is allowed to be larger: rrih < 160 GeV. 



In Figs. |3| |4] and [5| we have plotted the max(o"gj, Ugj) vs. for neutralino scattering with 
various restrictions. Note that these plots are made from independent scans and we have 
taken the scalar superpartners to be C(2 TeV). 

As discussed in Section [ll| we define "large" cross sections to be cig^^'^ > 5 x 10^^ pb and 
(TgQ^° > 10~^ pb, motivated by the projected near term range of current DD experiments. 
Hence, the shaded region delineates the (very approximate) reach of the next generation of 
SI and SD experiments. Note that this neglects the dependence of the sensitivity on the 



mass of the DM. The maximum for ctsd in Fig- 3 is given by Eq. (31) and for Figs. 4 and 5 



is given by Eq. (35). 



In Fig. [3] we show points for both thermal and non-thermal neutralinos. This is our most 
general framework, and in this case it is clear that the correlation between the relevant mixing 
angles (and hence cross sections) is weak. By only allowing points which have a thermal 
relic density within ±3 a of the WMAP measurement (see Figs. |4]and|5]), the correlation 
progressively improves. We will discuss this in detail in what follows. 

We will pay special attention to the tuh 00 limit. In any theory with multiple Higgs 
bosons, a small SI signal can occur when the diagrams from Higgs boson exchange cancel 
against one another. Two important points should be made. First, this cancellation is 
often incomplete and typically cannot be realized for scattering off of both protons and 
neutrons simultaneously. Second, such a cancellation is a conspiracy - it requires unexpected 
relationships between parameters in the Higgs sector and nuclear matrix elements. The finer 
the cancellation, the greater the conspiracy (for further discussion of this cancellation, see 



Sec. VIC). If one takes the decoupling limit for Fig. 3 so that SI DD is determined by 



h exchange alone, the maximum SI cross section is ~ 3 x 10 ^ pb. Note that even for 



niA ~ 0{TeV) there can be nontrivial contributions for tjs ~ (9(50) (see Eq. (22)). 



There is a negative correlation between fine-tuning and the size of DD cross sections (see 
Eq. (|26|) [Sni [60]. To emphasize this point, in Figs. [3| |4]and|5]we have marked points with 
< 500 GeV by blue dots and points with |/i| > 500 GeV by red crosses. The apparent 



feature around asi ~ 10~ pb in Fig. |4|is due to the finite range of rriA taken in this scan 
{ttla < 1 TeV) - the points above this gap have constructive contributions from h and H 
while the points below have destructive contributions. There are a few interesting features 
in Fig. |5| The gap which extends along the entire plotted range of SD cross sections is due 
to a slight cancellation between the various contributions from the light Higgs boson (see 



Eqs. (19) and (20)) which can occur at finite tfs {tp < 50 in this scan). The small number 
of points around ctsd = 3 x 10~^pb is due to the opening of the top threshold (see Fig. [2]). 
The behavior around ctsd = 2 x 10"^ pb is due to the cross over from dominantly Bino to 
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dominantly Higgsino DM, which occurs around = 500 GeV. 



Relic Density Not Imposed 
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FIG. 3: The max(crgj, CTgj) vs. figj-, cross sections in pb for the MSSM. The dots (in blue) and 
crosses (in red) correspond to < 500 GeV and |^| > 500 GeV respectively. The horizontal 
(vertical) line refers to the projected sensitivity for the next generation of SI (SD) experiments. We 
have shaded the near-term probeable region. Note that we are neglecting the dependence of this 
sensitivity on the neutralino mass. We have not imposed the thermal relic density constraint - all 
points are taken to have /Odm = 0.3 GeV/cm^, regardless of thermal abundance. All sfermions have 
masses of 0(2 TeV). If one takes the decoupling limit, there is a maximum value for usd = 3 x 10^^ 
pb. 
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FIG. 4: The max(cjgj, cjgj) vs. fjgj-, cross sections in pb for the MSSM. We have imposed that 
the thermal abundance of the neutrahnos is within it 3 cr of the WMAP measurement. The dots 
(in blue) and crosses (in red) correspond to \iJ,\ < 500 GeV and \fi\ > 500 GeV respectively. 
The horizontal (vertical) line refers to the projected sensitivity for the next generation of SI (SD) 
experiments. We have shaded the near-term probeable region. Note that we are neglecting the 
dependence of this sensitivity on the neutralino mass. All sfermions have masses of 0(2 TeV). 



A. Large SI and Large SD 



To have non-zero SI and SD signals, a Bino-Higgsino, Wino-Higgsino or Bino-Wino-Higgsino 
mix is required. In fact, appreciable SI and large SD signals can be generated as long as 
the Higgsino fraction is larger than C(10%). Note that the \fj,\ < 500 GeV points, which 
correspond to less fine-tuning in mz, imply large SD signals. When the gaugino fraction is 
dominated by Wino rather than Bino, the relative size of g and g' gives a slight enhancement 
in the SI cross section. There can be further enhancement of the SI cross section if sgn(ZB) ^ 
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FIG. 5: The max(o"gj, Cgj) vs. dgj-, cross sections in pb for the MSSM with gaugino mass unification. 
We have imposed that the thermal abundance of the neutralinos is within it 3 cj of the WMAP 
measurement. We have taken the decoupling limit (m^ = 4TeV). The dots (in blue) and crosses 
(in red) correspond to < 500 GeV and |^| > 500 GeV respectively (see the text for a discussion). 
The horizontal (vertical) line refers to the projected sensitivity for the next generation of SI (SD) 
experiments. We have shaded the near-term probeable region. Note that we are neglecting the 
dependence of this sensitivity on the neutralino mass. All sfermions have masses of 0{2 TeV). 



sgn{Zw) (see Eq. (19)) which accounts for points with the largest SI values in Figs, [sjand 
|4j This cannot occur in models with unified gaugino masses, where M2 ~ 2 Mi. 

Large SI and SD signals occur as long as there is non-trivial gaugino content in the WIMP. 
Imposition of the thermal relic density constraint for > mw, ensures a minimum required 
Bino component. If one imposes the large SI and SD conditions, < 0.7 and < 0.85 

below and above the top threshold respectively. Note that the large SD requirement implies 
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that < 200 GeV (see Fig. [2]). Hence, the assumption of a thermal history is necessary to 
conclude that the neutralino is a Bino-Higgsino admixture, rather than Wino-Higgsino. 

In the next three subsections we will attempt to elucidate the difficulties one encounters 
when trying to suppress SI and/or SD. This will allow us to argue that large SI and SD DD 
signals are the generic prediction for a well-tempered MSSM neutralino, since suppression of 
either SI or SD or both requires doing some gymnastics. While future data may force these 
contortions upon us, we conjecture that if the DM is a well-tempered neutralino, it is likely 
to be discovered in the next generation of DD experiments. 



B. Small SI and Small SD 

There are two ways to suppress both SI and SD. The first is to make \Zh^\ = \ZhJ = 0, 
which is equivalent to the yU — ?■ oo limit. This limit leads to fine-tuning of the electroweak 
scale. To achieve the proper thermal relic abundance in this case requires a Bino-Wino mix. 
Note that the Bino and Wino only mix indirectly through the Higgsino. Therefore, two 
insertions of the mixing factor are required, and the resulting mixing is of size (m^//i)^. One 
can see the effects of this limit by inspecting the red crosses in Figs. [3} |4] and |5} The upper 
bound in Figs. [3] and |4] are from points which are either Bino/Higgsino or fully mixed states 
while the points with the smallest values for SI are due to either Wino/Bino neutralinos or 
the cancellations discussed in Sec. IVI Dl 

The second option is to take Mi 2 ^ yU. This will imply that Zb,w = 0, thereby suppressing 
SI DD, and \Zh^\ = \ZhJ so that SD DD is also zero. Reproducing the measured relic 
density then requires ~ 1 TeV. When one does impose the thermal relic density as a prior. 
Fig. [2] shows that for DM masses of (9(TeV), i.e. the region of dominantly Higgsino DM, 
the SD cross section ranges from (9(10^^ pb) to 0. Fig. [s] shows the corresponding SI cross 
sections for this range. The trend of SI and SD going to zero in this plot is due to the limit 
Ml 2 — 00. Thermal dark matter in either of these two limits {fi or Mi_2 — > 00) will have a 
finely-tuned electroweak scale. Note that for either pure Wino or pure Higgsino DM there 
is a 1-loop diagram which leads to an SI DD cross section of (9(10^^^ pb) or (9(10^^^ pb) 
and an SD DD cross section of (9(10~^pb) or O(10~^'^pb) for the Wino or Higgsino case 
respectively [M]. We neglect this tiny contribution in our numerical scans. 
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C. Large SI and Small SD 



There are points which have large SI and SD with a nearly maximal gaugino fraction. If 
one relaxes the requirement of large SD, then the gaugino fraction can be pushed to nearly 
100% while keeping the product Zb,w Zh^ ^ approximately fixed, which in turn keeps the SI 
cross section constant. The relic density constraint can still be satisfied since both Winos 
and Higgsinos annihilate to bosons with approximately the same rate. 

There is another way to have small SD while allowing large SI. In the context of the SDM, one 
can take A = A', i.e. t/j = 1 in the MSSM. From the SDM mass matrix (see Appendix [B|), one 
can see that mixing between 5* and D_ will vanish. Since the SD cross section is proportional 
to this mixing factor, Z£)_, it will be zero as well. This effect accounts for the empty region 
in Figs. |3] and |4] since we restricted tp > 5 in our numerical scans. 

For tan /3 ^ 1.5, we find that for o"si ~ 5 x 10"^ pb the smallest cross section for SD is 
(TsD ~ 10~^ pb. If one allows cxsi < 5 x 10~^ pb, then as {Zh^^I ~^ 0' crso/fsi — \Zh^J'^ — >■ 0. 
Hence, SD falls off faster than SI. However, this is the /x — oo limit which leads to fine-tuning 
as described above. 



D. Small SI and Large SD 

Large SD requires a well-tempered neutralino, which naively also leads to large SI DD. In 
this section we will enumerate the various options one has for suppressing SI signals. We 
will argue that all options require fine-tuning or numerical coincidences^. 

Here are the options for minimizing asi: 

1. One can make and rriH heavy; however rrih ~ 115 GeV in the MSSM in the absence 
of large fine-tunings. Even in SplitSUSY, rrih^ 160 GeV. 

2. Since Cu,d ~ {Zw — trnZs), i-e the Higgs couples to the Zino, one could attempt to 
restrict the DM to only be a photino-Higgsino admixture. In Appendix [C} we show 
that this is impossible when one restricts M2 by the LEP bound. 

3. One can tune (/^^^ + 2 ^ fiS) against (/jj^ + f^rV + AV) ^ by tuning the 
contribution from H against that from h. As we will discuss below, it is not possible 



Another possibility is that both SI and SD from exchange of the Z° and Higgs boson respectively are 
small. If there exist light squarks, they can give rise to large SD signals [17]. Cross section estimates from 
light squark exchange are discussed more in Appendix lAl 



22 



to precisely tune this quantity to zero simultaneously for the proton and the neutron 
(see Fig. [6]). However, an approximate realization of this condition is possible - this 
is the tuning that underlies large SD/small SI points in Figs. |3]and|4]and reported in 
the literature {e.g. [62]). 

4. One can tune the contribution from the proton against the contribution from the 
neutron. The cancellation would only hold for a specific element. Since all experiments 
do not use the same elements, we will not pursue this case further. 

In what follows, we minimize the SI cross section by tuning the contributions from the h 



and H against each other (point 3 above). From Eq. (20), this cancellation requires (in 
the decoupling/large tp limit) sgsi{ZH^) = sgn(ZH^). This condition for cancellations to be 
possible was first noted in [63]. Using DarkSUSY we have confirmed that this is a necessary 
condition, not just in this limit, but for any values of the pseudo-scalar Higgs mass (m^) 
and tp. This condition only occurs for certain signs of Mi, M2 and /x. If large SD/small 
SI were observed for neutralino DM, this would constrain the signs in the neutralino mass 
matrix. 

Let us estimate the maximum allowed suppression. To good approximation^, the best one 
can do is to tune away the coupling to (for example) the proton: 

^ = ^ ( ■^Td' + frl + §f /tg ^ ^ = _ fip) ~ _x 64 — (36) 



In order for Eq. (36) to have a guaranteed solution requires independent control of a. and 
mil- Since there is a non-trivial relationship between a and mu (both are determined by 
m^), our lower bound provides a conservative estimate. Using Eq. (|9]) to estimate Cg and 



plugging in the relationship between and q from Eq. (36) gives cigj = and 



„ 4 2 (A - Z)2 2 2 1 



( + 2 I /S) f^l - (/S + /ff + I ) ' (37) 
^ 8 X 10'^^ pb ( ^^^^^ ^ (Iy)^ (SI with cancellations). (38) 

This gives an estimate for how small SI can be, absent taking some of Mi, M2, /i — ?■ 00. The 
effects of the current uncertainties on the hadronic matrix elements described in Sec. Ill Al 



From Fig. [g] the absolute minimum of the total SI cross section occurs between the region where the 
coupling to the proton and neutron vanish. Therefore, the following analytic estimate will be off by a 
factor of a few. 



23 



can change the amount of cancellation allowed (the coefficient in Eq. (36)), altering the lower 



bound in Eq. (38) by C(50%). 



In Fig. [6] we show the SI cross section on the proton, the neutron and both as a function of 
rriA for a 93 GeV neutralino with a thermal relic density of f^DM = 0.1, cigj-, = 9 x 10^^ 
pb and crg^ = 6 x lO""* pb. One can clearly see that both contributions to SI DD cannot 
both be canceled simultaneously. At the minimum, cTg}™ = 3 x 10~^^ pb for niA = 751 GeV. 
For a shift in of ~ 5%, the cross section becomes ~ 2 x 10^^° pb - a change of almost 
2 full orders of magnitude. This emphasizes the delicacy of the cancellation. Other than in 
the limited region where the cancellation occurs, the entire range is probeable by the next 
generation of SI experiments. 

Numerically, we find that for ctsd > 10^^ pb, the smallest cxsi can be is C(lO^^^pb) where 



the suppression beyond the value in Eq. (38) is due to small mixing angles. 



Finally , we note that while these kinds of conspiracies are allowed, there is no reason to 
expect that the SUSY breaking parameters have anything to do with the nuclear matrix 
elements. We take this as evidence that such cancellations are unlikely. 



VII. CONCLUSIONS 



In this work we have explored the physics of SD DD with an emphasis on the correlations 
with SI experiments. In the process, we have determined some expectations for the SD 
cross sections. In particular, in the MSSM, (cTg^^"^) < 6 x 10^'^ pb without making any 
assumptions about the thermal history. Again, allowing for a non-trivial cosmic history, 
but imposing the unified gaugino mass condition, we find (cg^^"^) < 4 x 10^^ pb. Finally, 
('^SD^'^) < 2 X 10~^pb when a thermal relic density is imposed. These represent important 
targets for future experiments. If one includes the possibility of squark exchange, a SD cross 
section as high as 2 x 10~^ pb can be reached for a neutralino which has a thermal abundance 
by utilizing the squark pole [HllES]- We note that in the absence of light squarks, if SD cross 
sections larger than ~ 6 x 10~^ pb were observed, the DM would not be an MSSM neutralino. 
This would point to more exotic theories like the SDM or models with light mediators [M] . 
For models which reproduce the relic density, in the decoupling limit, and unified gaugino 
masses, a 1-ton COUPP-like experiment could probe the entire range of SD cross sections 
up to WIMP masses of O(lTeV). 

More generally, we have argued that given the experimental constraints from LEP, neutralino 
DM is likely to be well-tempered with possible signals for the next generation of SI and SD DD 
experiments. In fact, any model (such as the SDM) which interacts with the SM via a light 
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FIG. 6: Plot of the SI DD cross section for the neutralino scattering off of a proton (sohd), a 
neutron (dashed) and both (dotted) as a function of ruA- For reference, the size of the SD cross 
section is about 9 x 10~^ pb (proton) and 6 x 10~^ pb (neutron) and = 93 GeV. The thermal 
relic density is r^oM = 0.1. The minimum value for the total SI DD is cjgj™ = 3 x 10^^^ pb for 
niA = 751 GeV. By changing mA by 5%, the cross section becomes ~ 2 x 10^^*^ pb. For small niA 
the cross section is on the order of crgi ~ 10^'' pb and in the decoupling limit the cross section is on 
the order of fisi ~ 10~^ pb - the entire region where there are not any conspiratorial cancellations 
is within the reach of the next generation of SI experiments. 

Higgs boson can imply a signal in SI experiments and any model of Majorana fermions with 
non-trivial couplings to the can imply a signal in SD experiments. We have enumerated 
the ways to avoid these arguments. Since all of these options involve a numerical conspiracy 
or some new source of tuning, we take them to be disfavored. With available methods we 
should be able to probe the majority of the natural range for the SI and SD DD signals of 
both thermal and non-thermal neutralino DM. 
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Appendix A: Squark Contributions to Direct Detection 



The neutralino can scatter off of quarks via s-channel squark exchange, giving contributions 
to O^^ or O^^. Only squarks that couple to the light quarks (m, d, s) will be able to contribute 
to the SI and SD cross sections since only the light quarks have non-negligible nuclear matrix 
elements. 

A non-zero "left-right" squark mixture is required since SI scattering converts a left-handed 
quark into a right-handed quark. Though a Bino/Wino mixture maximizes the coupling 
between the quarks and the neutralino, the scattering cross section for a pure Bino is of the 
same order. 

If one makes the standard assumption that left-right squark mixing [i.e. a-terms) are 
proportional to Yukawa couplings, then the squark mixing angle is proportional to mq/ifiq. 
Therefore, all SI couplings will be proportional to a quark mass and there is no enhancement 
for the light squarks over Higgs boson exchange. The maximum cross section is 

(a^r'\xN^ XN)) = 6 X 10"^ pb f ^^^V, (A.l) 

\ / max \ Tflg I 

for a Bino-Wino mix. This is subdominant to the Higgs boson exchange contribution barring 



the cancellations discussed in Sec. VI D 



The maximum crgD(xP ~^ XV) contribution from squark exchange is for a "left-handed" 
up-type squark coupling to a pure Wino, due to the larger SU{2) gauge coupling: 

(4r\xP ^ XP)) = 3 X 10-^ pb ('^^O^y . (A.2) 

\ / max y nT'u J 

This is typically subdominant to the Z° contribution to SD DD. Thus, we will focus on the 
effects of exchange in our discussions of the expected SD cross section. 



^ If exceptionally large left-right in the squark sector is allowed (perhaps through abnormally large a-terms) 
a contribution to crsi(x-^ ~^ X^) of 0(Wr^ pb) may be obtained. 
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Appendix B: The Bino/Higgsino and Wino/Higgsino Limits 



In the hmit of large Mi (M2) the neutrahno is dominantly a Wino/Higgsino (Bino/Higgsino) 
admixture. We can explore this effective 3 state system using the SDM defined as (see 



Eq. (27) above): 



>CsDM ^ fioD D + \ hS D + \' h* S D + ^ 



(B.l) 



The resulting lightest eigenstate (x) is specified by 

X = ZsS + ZdD + Z5D. (B.2) 

Following [9j, it is useful to write this system in a basis defined by S and D± = ^{D± D). 
Note that the labels ± have nothing to do with electric charge. The mass matrix is then, in 
the {S, D+, D_) basis, 



/ 



M 



SDM 



/is 75(A + A')t; ^(A-A')t; 
^ (A + A') V fiD 
}=J\-\')v -flD 



\ 



\ V2 

with the resulting lightest eigenstate, 

X = ZsS + Zd^ D+ + Zd_D_. 



(B.3) 



Since we are interested in the SD DD cross section, our goal is to extract the coupling of x 
to the Z". The coefficient of the operator of Eq. (11 ) is given by 

d, = --^\2Zo^ZoSTl (B.4) 

Note that \2 Zr,^ Z]:,_\ = |Z_dP — I^dP- Oiis can find analytic expressions for the mass 
eigenstates and the combination |2Z£i^ Zd_| in various useful limits. To second order in v, 
for l/i^l, \fis\, {\^J'D\ - |/is|) > Af, X'v 



1 2 Zo: Zo_ 



and for \fi£)\ = \fis\ ^ Af, X' v, 



\2Z 



D+ ^D- I 



(Y^ - A^) v"" 

2 2 ' 



(A' -X)v (A'2 - A2) 



2V2\fXs\ 



(B.5) 
(B.6) 

(B.7) 
(B.8) 
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Perturbing away from the limit of exact degeneracy gives corrections to these expressions of 
0{{fis — fJ'D)/f^D)- Note we have assumed that there is no CP violation for simplicity. In 
order to apply these expressions to the MSSM one can make the identifications 



SDM 


Bino/Higgsino 


Wino/Higgsino 




Ml 


M2 








\v 


-mz Sw 




X'v 


-mz Sw sjs 


mz Cyj S/3 



where we neglect terms of 0(1/M2) for the Bino/Higgsino system and C(l/Mi) for the 
Wino/Higgsino system. 

Explicitly making the substitutions for the MSSM we have 



-"w I 



and 



\Zh, 



|2 1/7 |2 



C2/3 cj ml 



2V2|mI 



for 






|Mi 


> mz, M2 — )■ 00 


for 


IM2I, 




IM2I 


> mz, Ml 00 



„2 „2 \ „2 



2V2\n\ 



for 


\Mi\ 


= l/^l 


for 


IM2I 





(B.9) 



(B.IO) 



Appendix C: No-go Theorem for photino-Higgsino DM 



The neutralino mass matrix in the (7, Z, H^, Hy) basis is given by 

/ Mic^ + M2S^ {Ml - M2) Cy, Sy, -mzS2wCp mzS2wSi3 \ 



M 



{Ml - M2) Cyj Sy, Ml si + M2 cl mz C2w cp -mz C2w sp 

-mzS2wCi3 mzC2wC(3 -/i 

mzS2wSp -mzC2wSp -yU 



Is it possible to generate a large SD/Sl ratio by having DM which is only a mixture of photino 
and Higgsino? The Higgsino component is required for a non-trivial coupling to the and 
an admixture of photino (and not Zino) will allow — 7^ without introducing 
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a coupling to the Higgs. We show that current phenomenological bounds preclude this 
possibility. 

There are two potential options. The first is decoupling the Zino by making it heavy while 
tuning the photino mass to be ~ /i. This implies taking the limit where Mi and M2 are 
large while the combination Mi + M2 stays small, which requires sgn(Mi) 7^ sgn(M2). 
Then the Zino-photino mixing will go like (Mi — M2)/m^ > 0{1). Note that we are free to 
take Ml < mz to suppress this mixing, but due to the LEP bound on the chargino mass, 
M2 > mz- The second option is to try to eliminate the photino-Zino mixing by taking 
Ml = M2. Then the Zino and photino have the same mass and the Higgsino will mix 
with both, resulting in a DM state which is an equal admixture of all 4 gauge eigenstates. 
Therefore, a neutralino cannot be a mixture of only photino and Higgsino. 
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